Non-converging hysteresis cycles in random spin networks 
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Behavior of hysteretic trajectories for cyclical input is investigated as a function of the internal 
structure of a system modeled by the classical random network of binary spins. Different regimes 
of hysteretic behavior are discovered for different network connectivity and topology. Surprisingly, 
hysteretic trajectories which do not converge at all are observed. They are shown to be associated 
with the presence of specific topological elements in the network structure, particularly with the 
fully interconnected spin groups of size equal or greater than 4. 

PACS numbers: 75.60.Lr, 75.75.+a, 89.75.Fb 
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Although hysteresis is often illustrated by closed loops, 
there is no fundamental reason to expect that period- 
ically varying external parameters should produce any 
kind of closed hysteretic trajectories in general. In fact, 
hysteretic trajectories which do not converge are ob- 
served experimentally in many systems. It is known, for 
example, that repeated mechanical loading and unload- 
ing of some materials leads to gradual accumulation of 
strain and eventual material failure, rather than stable 
hysteresis cycles [![. Similar non-convergent behavior 
in magnetic materials is manifested as a gradual shift 
of magnetization occurring under the action of periodic 
external field [2j. The question arises: What in the 
structure of a system determines the existence of non- 
convergent hysteretic trajectories? 

Here this issue is investigated by looking at binary spin 
networks as a way to model the structure of a system. Al- 
though it is most frequently employed to model magnetic 
materials, binary spin network is a classical prototype for 
studying hysteretic properties of many natural systems. 
Hysteretic cycles have been studied previously using the 
Random Field Ising model Q and random networks of 
antiferromagnetically coupled spins Q where the focus 
was mainly on the effects of disorder. In addition, inter- 
esting memory effects have been discovered in spin glass 
networks such as complementary point and reversal-field 
memories Particularly important property found to 
occur in ferromagnetic [fj and some antiferromagnetic 
networks [7| is the return point memory (RPM), which is 
responsible for the recovery of the microscopic state and 
formation of closed minor hysteresis loops after the very 
first external field cycle. Networks without RPM have 
been found to display only a gradual minor loop forma- 
tion or sub-harmonic cycles Q (cycles with period multi- 
ple of the external field period). Despite the substantial 
amount of work performed in this area, ability of binary 
networks to describe non-convergent hysteretic trajecto- 
ries was not observed. The link between the network 
structure and the convergence character of hysteretic tra- 
jectories also remains largely unexplored. 

In this Letter, investigating the link between the net- 
work connectivity and topology and its hysteretic be- 
havior, we discover growth of non-convergent trajectories 



with emergence of certain topological structures. Specif- 
ically we find that RPM observed at low network con- 
nectivity gives way to transient trajectories converging 
to sub-harmonic cycles or to trajectories with very long 
transient time growing without bound as the system size 
increases and thus becoming non-convergent. 

As a prototypical system, we consider classical Erdos- 
Renyi (ER) random network of anti-ferromagnetically 
coupled spins. It is one the simplest realizations of a com- 
plex network. It is also convenient because its connectiv- 
ity and topological structure can be tuned by adjusting 
a single parameter - probability of connection between 
pairs of spins [9( . Generally, a random network of N spins 
can be viewed as a random graph with TV nodes (ver- 
tices) connected by undirected edges representing pair- 
wise interactions between the spins. ER networks are 
constructed using a random graph process where edges 
are added one by one between randomly chosen pairs of 
spins. The structure of the graph is characterized by 
its adjacency matrix Aij with elements Ay — 1 if the 
vertices i and j are connected or Ay = if the ver- 
tices are disconnected. Antiferromagnetic interactions 
between the spins will then be described by a matrix 
Jij = —J A^, where the J is a positive number repre- 
senting the interaction strength. 

Dynamics of the system is governed by minimization 
of the standard Ising-type energy: 
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where Sj denotes the ±1 spin variable and H is the exter- 
nal field. The first and the second sum in (Q]) are respec- 
tively the spin-spin and the spin-external field interaction 
energies. Additionally, there is a symmetric double-well 
potential energy /j, due to which any isolated spin flips 
from —1 to +1 at the field-threshold +£*,, while the flip- 
ping from +1 to —1 occurs at the field — <x;. The positive 
threshold magnitudes on will be viewed here as random 
variables to mimic a structural disorder in the system. 
Note that the difference between the 'up' and 'down' 
thresholds, +a, > —on, results in an inherent hysteretic 
behavior of each spin. Such classical hysteretic spins are 
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frequently used as representations of single-domain mag- 
netic grains, tiny capillary pores in absorbing materials, 
vortex pinning imperfections in superconductors, indi- 
vidual decision making agents in socio-economic systems, 
etc. They are also the main building blocks of some phe- 
nomenological models of hysteresis pjj] . 

The local field hi governing the flipping of each spin is 
obtained as a variation of Eq. ([T]) with respect to and 
equals: 



hi = -SG/Ssi = -jJ2 A 



H. 



(2) 



We consider standard field driven adiabatic dynamics 
[f| where, at each time step, the spin states are up- 
dated only if hiSi < —on. For a network with nega- 
tive interactions such dynamics may produce avalanches 
with backward flipping spins when AHAsi < result- 
ing in a non-monotonic variation of state of the entire 
system even when the field H varies monotonically 
Here we avoid such backward spin flips by restricting 
dynamic behavior to single spin flips. According to (|2|) 
the spin back-flips do not occur as long as the inequality 
J ■ max(di) < rnSn{ai) — a m in is satisfied for the partic- 
ular network, where di = J2j A ij is a degree of vertex i 
in the network, i.e. the coordination number of the spin 
Si. Hence, choosing J < J t — a m - m / max(di) guarantees 
the absence of back- flips, resulting in a simple avalanche 
behavior consisting of single spin transitions AHAsi > 
and monotonic state variation. In simulations, we use 
Gaussian distribution of ai with variance a and mean 
fi » a 12|. Such an assumption is natural for various 



realistic systems e.g. magnetic films with very stron g pe r- 
pendicular anisotropy or patterned nanostructures [13[ . 

It turns out that even the single spin-flip dynamics 
limit introduced above yields complex minor loop behav- 
ior. This is demonstrated here by simulating hysteretic 
trajectories corresponding to periodic cycles of the exter- 
nal field (Fig. 1) obtained by increasing the field starting 
from the negative saturation, where all spins are in the 
— 1 state, to the point H r with the average spin state 
(magnetization) (S), then decreasing the field to — H r 
and returning back to H r . The ability of the network 
to recover its state when the field returns to H r will be 
quantified by the cycle opening, denoted by AS and equal 
to %-difference in a number of spins which did not return 
to the original state. 

We observe that the cycle opening AS* depends on the 
degree d (i.e. connectivity) of the network, which is de- 
fined as an average over all di in the particular network 
realization. For connectivity d « 1, the system returns 
to the same state at the end of the very first cycle. This 
is expected since the majority of spins is isolated or form 
couples and such antiferromagnetic systems are known to 
have a RPM [ljj]. As d increases the cycle opening AS 
becomes nonzero. However, a well pronounced increase of 
AS is observed only after the percolation threshold of the 
network at d « 1, when a giant spin cluster appears Q. 
Then, after reaching a maximum, AS starts to decrease 
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FIG. 1: Hysteretic cycles for two periods of external field. S 
and (5*) denote, respectively, the spin microstate state and 
the average state (magnetization). AS is a measure of the 
cycle opening calculated by comparing spin patterns before 
and after each cycle. Results were obtained from simulation 
of 50 2 random spin-network with the average coordination 
number (degree) d = (di) — 10. Reducing the degree to 
d < 1 results in the immediately closed minor loops (AS = 0) . 
Inset: Dependence of AS on the external field period number. 
The steady state reached after initial transient time tr can 
contain simple minor loops with ts = 1 or subharmonic cycles 
with ts > 1. 



to zero as d approaches the network size N with the limit 
AS = obtained for fully connected network d = N. 
Such behavior is expected since the fully connected spin 
network can be viewed as a mean-field Preisach model, 
which has been shown to have RPM [To] ]. 

Numerical tests for different network sizes N (up to 
50 2 ), different disorder a and different average energy 
per spin A = (JYlj A ij s j) ~ J ■ d - (S), revealed that 
the cycle opening AS" depends on the ratio Aja rather 
than on A and a separately. In addition, AS was ob- 
served to be independent of the network size as long as 
d << N, and we found no dependence on the mean /j, 
in the assumed limit /i >> a. Since the probability of 
finding various topological interconnection structures in 
the ER network, such as trees and cliques (completely 
interconnected sub-graphs) @, depends both on d and 
N, these results demonstrate that the first cycle open- 
ing does not depend on the topological properties of the 
network. We note, that this conclusion is supported by 
similar results obtained for a regular lattice binary spin 
models, where AS was calculated as a function of the in- 
teraction range (which is. of course, proportional to the 
coordination number) [161 ]. 

It turns out, however, that the network topology does 
determine the minor loop formation behavior. To show 
this we studied the dependence of AS* on the number 
of external field cycles, the network size N, the Aja ra- 
tio and the connectivity parameter d. An example of 
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FIG. 2: a)-d) contour maps showing the cycle opening AS 
for different values of interaction energy A/ct and the network 
degree d (note the logarithmic scale of A/a and d axes) for 
respectively: 1-st, 10-th, 50-th and 100-th cycle. The lower 
bounds for the 'limiting' region with non-converging loops 
correspond to (A/a)t ~ 2.3 and d t ~ 13. All minor loops have 
been obtained by reversing the field at H r where (S) = 0.2, 
when the effect was most pronounced. Data are for N = 50 2 , 
ct = 0.1 and averaged over 50 random graph and disorder 
realizations. Error bars level is about 1%. 



the AS (A/ a, d) function is given by the contour plots in 
the Fig. 2 a-d for four subsequent field cycle numbers 
and a fixed network size N = 50 2 . The figure illustrates 
that in the low (A/ a, ci)-parameter region closed minor 
loops with AS — appear already after a few initial field 
cycles (Fig. 2 a-b). Similar behavior, termed as 'tilting 
effect' (bascule), is frequently observed in magnetism [2j 
especially for clean ferromagnets. It is often attributed to 
the coupling between only a small numbers of magnetic 
domains (consistently with the low connectivity network 
picture considered here). 

For the parameter region bounded from below by cer- 
tain critical values (A/er) t and d t (Fig. 2 d) the behavior 
changes dramatically and closed minor loops often do not 
form even after 100 field periods. Quite surprisingly, we 
find that the value dt is remarkably close to the theoret- 
ical value at which the ER network is known to undergo 
a topological transition associated with the emergence of 
cliques of size 4 (fully interconnected groups of 4 spins) 
Q . Two new types of cyclic behaviors are observed: cy- 
cles with a very long transient time, tx > 100, and sub- 
harmonic cycles with a spin pattern recovering repeatedly 
with period ts > 1 (tr an d ts are defined in the Fig. 1). 
Note that subharmonic cycles were also observed previ- 
ously for some spin-glass networks p|. Here, the period 
ts is found to grow with the number of 4-cliques present 
in the network. However, this growth is substantially 
slower than the increase of the network size N and no 
conclusions can be made regarding the behavior of the 



FIG. 3: a) Dependence of the transient length on the number 
of 4-cliques for N = 20 2 and 10 3 . (tr) and (C4) are averages 
over 100 network realizations, b) {tr) as a function of N ob- 
tained for (C4) = 0.03, 0.3 and 2. Lines for different (C4) are 
exponential fits. Data correspond to reversal magnetization 
(S) = 0.2, A/ct = 10, and ct = 0.1. Inset: dependence of the 
exponential fit parameter r on (C4). r grows without bound 
with diminishing number of 4-cliques C4. 



sub-harmonic cycle length in the thermodynamic limit 
(JV->oo). 

On the other hand, the transient length grows rapidly 
as shown in the Fig. 3a by the log-log dependence of 
average (tr) on the average number of 4-cliques, (C4), 
(averaged over 100 network realizations) for sufficiently 
strong interaction. Dependence on the network size is 
plotted in the Fig. 3b and follows the exponential law 
(tr) oc exp(iV/r). The parameter r has been obtained 
by fitting and, according to the inset in Fig. 3c, r — > 00 
in the limit (C4) — * 0. This means that for d > dt, when 
4-cliques emerge, the (tr) grows faster than TV, implying 
that hysteretic trajectories will not converge in the ther- 
modynamic limit (N -> 00). 

It is surprising that the present antiferromagnetic bi- 
nary spin network model with a simple single spin-flip 
field driven dynamics displays several qualitatively very 
different regimes, including a regime where hysteretic tra- 
jectories do not converge. Assuming that interactions 
are strong enough compared with the disorder, the de- 
termining factors appear to be the network connectivity 
and its topology. This suggests possible design rules for 
obtaining desired hysteretic behavior in artificial mate- 
rials. Moreover, this conclusion suggests that behavior 
of hysteretic cycles could prove to be a useful charac- 
terization method for probing topology and connectivity 
of some systems with complex interactions. Care is re- 
quired in carrying out such characterization because, due 
to noise, non-convergent cyclic behavior can be easily 
missed if only macroscopic measurements are performed. 
For example, average state (magnetization) in our sim- 
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ulations did not show the same behavior as the micro- 
scopic state. It is also interesting to note, that the model 
employed here did not reproduce some other forms of 
non-convergent behavior such as the reptation effect 
observed in some magnetic systems. It is possible, that 
the reptation effects will be observable only in spin net- 



works with mixed, i.e. ferromagnetic and antiferromag- 
netic interactions or in networks where spins have vector 
nature. 
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